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Measurements of the total and differential cross sections of Higgs boson production are performed 
using 20.3 fb~^ of pp collisions produced by the Large Hadron Collider at a center-of-mass energy 
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H > 77 and H —>■ ZZ* —>■ event yields, which are combined accounting for detector efficiencies, 
fiducial acceptances and branching fractions. Differential cross sections are reported as a function of 
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33.0 ± 5.3 (stat) ± 1.6 (sys) pb. The measurements are compared to state-of-the-art predictions. 
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This Letter presents measurements of the total and 
differential cross sections of inclusive Higgs boson pro¬ 
duction using 20.3 fb“^ of pp collisions produced by the 
Large Hadron Collider (LHC) [1] at a center-of-mass en¬ 
ergy of y/s = 8 TeV and recorded by the ATLAS detec¬ 
tor [2] . The measured cross sections probe the properties 
of the Higgs boson and can be directly compared to the 
theoretical modeling of different Higgs boson production 
mechanisms, such as the most recent gluon fusion (ggF) 
QCD calculations. They can also be used to constrain 
new physics scenarios, for example using the effective 
field theory framework as proposed in Refs. [3-7]. The 
analysis uses event yields measured in the iL —>■ 77 and 
iL —>■ ZZ* —>• 4£ decays and detector efficiencies, both de¬ 
termined as described in Refs. [ 8 , 9]. The statistical un¬ 
certainties on the Higgs boson signal yields in both chan¬ 
nels are larger than the systematic uncertainties, while 
the total uncertainties in the two channels are similar. 
Combining the analyses improves the precision of the 
cross-section measurements by up to 40%, and by 25-30% 
on average, with respect to the corresponding measure¬ 
ments in the most precise individual channel. 

Distributions of the differential pp ^ H cross sections 
are reported as a function of the transverse momentum 
p^ and the rapidity jj/^j of the Higgs boson, the jet multi¬ 
plicity A^jets, and the transverse momentum of the lead¬ 
ing jet pJjl. The observables p^ and jy^j describe the 
kinematics of the Higgs boson. They are sensitive to 
perturbative QCD modeling in ggF production, which is 
the dominant Higgs boson production mechanism in the 
Standard Model (SM). The jy^j distribution furthermore 
offers a clean probe of the gluon parton distribution func¬ 
tion (PDF) and will play a role in future PDF fits. The 
Ajets and observables probe the theoretical modeling 
of partonic radiation in ggF production as well as the 
overall rate and modeling of jets in vector-boson fusion 
(VBF) and associated Higgs boson production {VH and 
tiH). Jets produced in VBF, VH and tiH processes tend 


to have higher transverse momenta than those produced 
via ggF production, however the sensitivity to measuring 
these contributions is weak with the current amount of 
data. 

Cross sections are extracted using a combined like¬ 
lihood built from the signal yields in the iL —> 77 
channel and the data and background yields in the 
H —>■ ZZ* —>■ 4£ channel, as well as detector efficiencies, 
fiducial acceptances and SM branching fractions [10]. A 
complementary approach, using a separate likelihood, 
measures the shape of the differential distributions by 
imposing a unity normalization constraint, which re¬ 
moves the implicit SM assumption on the branching frac¬ 
tions. For the extraction of the signal yields and the cor¬ 
rections of detector efficiencies, it is assumed that the 
signal in both channels is due to a narrow resonance 
with a mass mn = 125.36 ± 0.41 GeV as measured by 
the ATLAS Collaboration [11]. The signal yield in the 
iL —>■ 77 channel is obtained from fits to the diphoton 
mass spectra [ 8 ], and from the background subtracted 
data yield in a mu mass window of 118 to 129 GeV for 
the H ^ ZZ* ^ M channel [9]. The fiducial acceptance 
in both channels [ 8 , 9] is derived using a set of Monte 
Carlo (MC) event generators. Powheg-box [12-14], in¬ 
terfaced with Pythia8 [15] for showering, is used to gen¬ 
erate ggF and VBF events, while Pythia8 is used to 
simulate VH and associated production with top quarks 
(ttH) and 6 -quarks (bbH). The fiducial acceptance for 
events with |y^| < 1.2 is approximately 72% for H —)• 77 , 
and 55-59% for JL —>■ ZZ* —>■ For higher |y^|, the ac¬ 
ceptance decreases to 35-38% in both channels. The fidu¬ 
cial acceptance is more constant as a function of the other 
variables and is in the range 56-62% for the iL —?► 77 
channel and 44-53% for the H ^ ZZ* ^ channel. 

After correcting the differential cross sections and nor¬ 
malized shapes for fiducial acceptance and branching 
fractions, the corresponding measurements in both chan¬ 
nels are found to be in good agreement with each other; 
p-values obtained from compatibility tests are in the 
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range 56-99%. 

In the binned maximum-likelihood fit, the statisti¬ 
cal uncertainty of the H —^ 77 event yield is modeled 
using a Gaussian distribution, while the event yield 
in the H ^ ZZ* ^ 4£ channel follows a Poisson dis¬ 
tribution due to the small sample size. Experimen¬ 
tal and theoretical systematic uncertainties affecting the 
signal yields, detector efficiencies, branching fractions 
and fiducial acceptance corrections are taken into ac¬ 
count in the likelihood as constrained nuisance param¬ 
eters. Nuisance parameters describing the same uncer¬ 
tainty sources are treated as fully correlated between 
bins and channels. Systematic uncertainties on the 
i/ —>■ 77 and H ^ ZZ* ^ ii background estimates and 
efficiency correction factors, as well as the uncertainty 
on the integrated luminosity, are described in detail in 
Refs. [ 8 , 9]. The branching fraction uncertainty due to 
the assumed quark masses and other theoretical uncer¬ 
tainties are evaluated following the recommendations of 
Ref. [16], considering uncertainty correlations between 
the i/ —>■ 77 and H —>■ ZZ* —>■ A£ decay channels. Un¬ 
certainties on the acceptance correction related to the 
choice of PDF set are evaluated by taking the envelope 
of the sum in quadratures of eigenvector variations of 
the baseline (CTIO [17]) and the central values of alter¬ 
native (MSTW2008NLO [18] and NNPDF2.3 [19]) PDF 
sets. Uncertainties on the acceptance correction asso¬ 
ciated with missing higher-order corrections are evalu¬ 
ated by varying the renormalization and factorization 
scales coherently and individually by factors of 0.5 and 
2 from their nominal values, and by reweighting the 
distribution from Powheg-box to the prediction of the 
HRes 2.2 calculation [20, 21]. The envelope of the max¬ 
imum deviation of the combined scale variations and the 
reweighting is used as the systematic variation. To 
account for the uncertainty in the mass measurement, 
the Higgs boson mass is varied by ±0.4 GeV. To as¬ 
sess the systematic uncertainty due to the assumption of 
SM cross-section fractions of the Higgs boson production 
modes, the VBF and VH fractions are varied by factors of 
0.5 and 2 from the SM prediction and the fraction of tiH 
is varied by factors of 0 and 5. These factors are based 
on current experimental bounds [22-26]. The total un¬ 
certainties on the acceptance correction range from 1 % 
to 6 %, depending on the channel, distribution and bin . 

The total systematic uncertainties on the combined dif¬ 
ferential cross sections range from 4% to 12%, depending 
on the distribution and bin. For the kinematic variables 
Pt and Ip^I, the largest systematic uncertainties on the 
differential cross sections are due to the luminosity and 
the background estimates in both channels. For the jet 
variables fVjets and pi ^, the largest systematic uncertain¬ 
ties on the differential cross sections are due to the jet en¬ 
ergy scale and resolution. In the shape combination, the 
normalization uncertainties including luminosity, branch¬ 
ing fractions, and efficiency uncertainties do not apply. 
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FIG. 1. Measured total cross section of Higgs boson produc¬ 
tion compared to two calculations of the ggF cross section. 
Contributions from other relevant Higgs boson production 
modes (VBF, VH, tiH, bbH) are added using cross sections 
and uncertainties from Ref. [10]. Details of the predictions 
are presented in Table I. 


ATLAS \ls = 8TeV, 20.3 fb' 

^ ^ /-/^ZZ*->4/ 

{ comb, data ■ syst. unc. 




pp-^H, 125.4 GeV 

+ a;,H = 3'0±0-1 Pb 
XH = VBF +VH+ tlH + bbH 

■ QCD scale uncertainty 

■ Total uncertainty (scale © pdf+q^) 


NNLO+NNLL 

LHC-XS 


N^LO 

ADDFGHLM 


Statistical uncertainties dominate all resulting distribu¬ 
tions, ranging from 23% to 75%. 

The total pp ^ H cross section is determined in the 


TABLE I. Summary of the ggF predictions used in the 
comparison with the measured cross sections. The second 
column states the order in QCD perturbation theory and 
which threshold resummation is applied, if any. Further de¬ 
tails are provided in the footnotes. All predictions are for 
tuh = 125.4 GeV and Ys = 8 TeV. 


Total cross-section calculations 

LHC-XS [10] 
ADDFGHLM [27-30] 

NNLO±NNLL 

N3lQ 

Analytical differential 

cross-section predictions 

HRes 2.2 [20, 21] 

STWZ [31], BLPTW [32] 
JetVHeto 2.0 [33-35] 

NNLO±NNLL 

NNLOTNNLL 

NNLOTNNLL 

Monte Carlo event generators 

SHERPA 2.1.1 [36, 37] 
MG5_aMC@NLO [38, 39] 
POWHEG NnLOPS [40, 41] 

H± 0,1, 2 jets @NLO ‘’j 
i7±0, l,2jets @NLO 
NNLO>oi, NLO>'=i‘j”' 


Considers b- (and c-) quark masses in the gg H loop 
^ Includes electroweak corrections 
^ Based on MSTW2008nnlo [18] (oa from PDF set) 

Uses TT^ -resummed gg ^ H form factor 
® NNLO refers to the total cross section 
^ Based on the CTlOnnlo PDF set 

s In the notation of Ref. [31], this corresponds to NNLL' 

^ Includes 1-jet resummation included at NLL'H-NLO 
^ Based on the CTlOnlo PDF set 

j Uses MEPS@NLO method and CKKW merging scheme [42-44] 
^ Software version 2.2.1, NLO merged using FxFx scheme [39] 

^ Interfaced with PythiaS for parton showering 
^ Uses Minlo method & reweighting to HNNLO [41, 45, 46]. 
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i/ — >■ 77 channel to be 31.4 ± 7.2 (stat) ± 1.6 (sys) pb and 
in the H ZZ* —>■ 4£ channel to be 35.0 ± 8.4 (stat) ± 
1.8 (sys) pb. Combining the analyses yields Upp^H = 
33.0 ± 5.3 (stat) ± 1.6 (sys) pb. Figure 1 presents a com¬ 
parison of these measurements with two ggF predictions 
to which contributions from other relevant Higgs boson 
production modes (VBF, VH, ttH, hbH) are added us¬ 
ing cross sections and uncertainties from Ref. [10]. The 
LHC-XS ggF prediction, recommended in Ref. [10], is ac¬ 
curate to next-to-next-to-leading order (NNLO) in QCD 
and utilises threshold resummation accurate to next-to- 
next-to-leading logarithms (NNLL). A signihcant effort 
has been undertaken by the theory community to provide 
ggF cross sections beyond this precision through various 
improvements in the perturbative calculations [31, 47- 
51]. Recently, the ADDFGHLM group has provided a 
fixed-order calculation accurate to next-to-next-to-next- 
leading order (N^LO) [27-30]. A PDF uncertainty of 
^g;g% is assigned to the LHC-XS prediction, derived fol¬ 
lowing the recommendations in Ref. [16]. This uncer¬ 
tainty is increased by 11q'i% for the ADDFGHLM pre¬ 
diction corresponding to the change in uncertainty of the 
MSTW2008nnlo PDF set when changing the calculation 
from NNLO to N^LO. The PDF uncertainty is treated 
as uncorrelated with the QCD scale uncertainty. 

The central value of the measured total cross section 
is larger than the SM predictions presented in Fig. 1. 
A likelihood-ratio test statistic is used to quantify the 
agreement, using a bifurcated Gaussian to model the 
asymmetric theory uncertainties. The resulting p-values 
are 5.5% and 9.0% for the agreement between data and 
the predictions from LHC-XS and ADDFGHLM, respec¬ 
tively. The ratio of the measured cross section to the 
LHC-XS prediction is higher than the results presented 
in Refs. [22, 23, 58], which use an event categorization 
based on the expected SM yields in the different Higgs 
boson production modes. 

The larger Higgs event yield observed in data moti¬ 
vates measurements of differential cross sections to in¬ 
vestigate if the excess is localized to specific kinematic 
regions. Figure 2 shows the comparison of the combined 
cross sections in different inclusive and exclusive jet mul¬ 
tiplicity bins with state-of-the-art predictions, including 
NLO-accurate multi-leg (ML) merged ggF MC event gen¬ 
erators (further details are given in Table I). Jets are 
reconstructed using the anti-Zc^ algorithm [52] with a ra¬ 
dius parameter R — 0.4 [53], and are required to have 
Pt > 30 GeV and \y\ < 4.4. Simulated particle-level 
jets are built from all particles with cr > 10 mm exclud¬ 
ing neutrinos, electrons and muons that do not originate 
from hadronic decays. Photons are excluded from jet¬ 
finding if they lie inside a cone of radius AR < 0.1 of 
an electron or muon, and neither the photon nor lep¬ 
ton originate from a hadron decay. To allow compar¬ 
isons with the unfolded measurements, the analytical cal¬ 
culations are corrected for effects of hadronization and 



FIG. 2. Measured Higgs boson production cross sections 
in inclusive and exclusive jet multiplicity bins compared to 
different theoretical predictions (see Table I for details and 
references). 


multiple particle interactions. These correction factors 
and their associated uncertainties are obtained using the 
Pythia8 and Herwig [54] MG event generators with 
different tunes [55-57] . The total cross sections from the 
ML merged predictions are lower than from fully inclusive 
NNLO-I-NNLL calculations. However, for Ajets > 1, the 
MC predictions formally have NLO accuracy, which is the 
same as the analytical calculations. Contributions from 
other relevant Higgs boson production modes are gener¬ 
ated using PowHEG for VBF and PythiaS for VH, ttH, 
and bbH, and are scaled to the cross sections in Ref. [10]. 
Uncertainties are assigned to all MC predictions from 
QCD scale and PDF variations. The ML-merged ggF 
predictions also have uncertainties due to the choice of 
merging scale. The SHERPA uncertainties further in¬ 
clude resummation scale variations. The measured cross 
sections are higher than the predictions for all measured 
jet multiplicities. The poorest agreement between data 
and predictions can be found in the inclusive and exclu¬ 
sive 1-jet bins, with local p-values ranging between 0.1% 
and 3.6%. Normalizing the total expected cross section 
to the data results in an improved agreement for these 
bins, with local p-values ranging from 4-29%. 

The combined differential cross sections as a function 
of p^, Jy^l, and are shown in Fig. 3 (left). The 
measured p^ and distributions are compared to 
the HRes calculation and the pi^ measurement is com¬ 
pared to STWZ and JetVHeto predictions. Figure 3 
(right) shows the comparisons of the normalized shapes 
to predictions from the MC event generators NNLOPS, 
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FIG. 3. Differential cross sections (left) and normalized cross-section shapes (right) for inclusive Higgs boson production 
measured by combining the H —>■ 77 and H —>■ ZZ* —>■ 4£ channels. The measured variables are the Higgs boson transverse 
momentum pJj! (top) and its rapidity |j/^| (middle), and the transverse momentum of the leading jet (bottom). The 0-30 GeV 
bin of the distributions corresponds to events without jets above 30 GeV. Various theoretical predictions are presented, 
using the same bin widths as the measurement. 
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SHERPA 2.1.1, and MG5_aMC@NL0, as well as the 
HRes calculation. The uncertainties on the predicted 
shapes are evaluated following the same approach as for 
the differential cross-section predictions. They are de¬ 
rived from the impact of QCD scale, merging scale and 
PDF variations. The mean of the measured distri¬ 
bution is 40.1 ± 3.0 GeV, while the means of the MG 
predictions range from 34 to 37 GeV. 

The p-values quantifying the compatibility of the mea¬ 
sured cross sections and predictions range from 2 % to 
26%, and for the shapes from 8 % to 88 %. For the calcula¬ 
tion of these values, the theory uncertainties are assumed 
to be Gaussian distributed and fully correlated between 
bins. 

In conclusion, this Letter presents the first measure¬ 
ments of total and differential cross sections and shapes 
for inclusive pp ^ H production. The measurements 
were performed in the iL —?► 77 and H —>■ ZZ* —>■ 4£ 
channels using the full 2012 dataset, which consists of 
20.3 fb“^ of pp collisions produced by the LHC at a 
center-of-mass energy of yi = 8 TeV and recorded by 
the ATLAS detector. The results of the two channels 
are compatible and have similar precision. The measure¬ 
ments indicate that the total production cross section of 
the Higgs boson is larger, and that it is produced with 
larger transverse momentum and more associated jets 
than predicted by the current most advanced SM cal¬ 
culations, however more data is needed to confirm these 
observations. 
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SUPPLEMENTAL MATERIAL 


Fiducial acceptance 


The fiducial cross section ai in a given bin i can be 
expressed as 

^2 f* Y 2 Y 2 ^ ^ ' 

L D at Ci o OLi 

where is the measured Higgs boson signal yield, £ is 
the integrated luminosity (20.3 fb“^ for this analysis), B 
is the branching ratio (0.228% for i? —)■ 77 and 0.0129% 
ioT H ^ ZZ* ^ I = e or /it), Oi is the fiducial accep¬ 
tance and Ci is a correction factor for detector effects, pri¬ 
marily accounting for reconstruction efficiency but also 
for bin-to-bin migration. For H ZZ* —>■ 4^, the signal 
yield is defined as the number of observed events Udata 
in a window around the Higgs boson mass peak minus 
the background estimate: rii = ndata.i — fibkg, 2 , while 
for H —>■ 77 , the signal yield is extracted from a simul¬ 
taneous signal-bbackground fit of the distribution. 
The correction factors for detector effects Ci, along with 
their systematic uncertainties are taken from the differ¬ 
ential cross section measurements in the individual chan¬ 
nels [ 8 , 9]. The differential cross section is defined as the 
fiducial cross section divided by the bin width. 


For each bin i, the acceptance factor for each decay 
channel is defined as 


^incl,2 

The fiducial acceptances for both channels and all mea¬ 
sured distributions are presented in Fig. 4 and Tables H 
and HI. They are based on Eq. 2 and derived using the 
Higgs MC samples described in the text. Forp^ and |y^|, 
ai is the probability for an event to pass the fiducial re¬ 
quirements. The acceptance is lower for H —>■ ZZ* —>■ M 
than for iJ —> 77 since it is less likely for four decay prod¬ 
ucts to fulfill the fiducial requirements. For the jet vari¬ 
ables and iVjets) an additional migration effect en¬ 
ters due to overlap between jets and the Higgs boson de¬ 
cay products, which affects the fiducial regions differently 
than the total phase space, where no Higgs boson decay 
products need to be considered. The fiducial acceptance 
falls off steeply as the Higgs boson rapidity increases, as 
both hducial definitions include pseudo-rapidity require¬ 
ments on the Higgs boson decay products. 

Additional figures 


TABLE II. Fiducial acceptance factors in percent for the 
H —> ZZ* —>■ M measurement with associated uncertainties. 
The binning is the same as in Fig. 4. 


Bin 

1 

2 

3 

4 

5 

Incl. 

-^jets 

46.7 ± 1.1 

45.0 ± 1.1 

47.8 ± 1.0 

49.8 ± 1.3 

50.0 ± 1.9 


ph 

45.0 ± 1.1 

46.2 ±0.8 

47.7 ± 0.8 

50.3 ± 0.7 


\y^\ 

59.2 ±0.9 

58.7 ± 0.7 

57.9 ±0.8 

55.3 ±0.6 

34.8 ± 0.8 

p? 

44.4 ± 1.0 

45.2 ±0.9 

47.6 ±0.7 

52.9 ± 0.7 



Figure 5 presents the measured jet multiplicity distri¬ 
butions. The lower two subfigures include the individual 
77 and H —)■ ZZ* —measurements. Figure 6 
presents the same six distributions as shown in Fig. 3, 
but with the individual channel measurements overlaid. 


TABLE III. Fiducial acceptance factors in percent for the 
R —> 77 measurement with associated uncertainties. The bin¬ 
ning is the same as in Fig. 4. 


Bin 

1 

2 

3 

4 

5 

6 

7 

8 

Incl. 

60.4 ± 1.2 








Ajets 

62.1 ± 1.4 

59.1 ±0.7 

58.4 ± 1.4 

55.9 ±3.5 





pipi 

62.1 ± 1.4 

58.9 ±0.9 

57.4 ±0.9 

57.2 ±0.9 

58.3 ± 1.6 




\y^\ 

71.5 ± 1.2 

71.6 ± 1.2 

71.7± 1.2 

71.5 ± 1.2 

67.7 ± 1.3 

37.7 ±0.5 



Pt 

62.2 ± 1.1 

62.2 ± 1.2 

60.6 ± 1.3 

59.3 ± 1.1 

57.9 ± 1.2 

57.5 ±0.9 

56.7 ± 1.1 

58.6 ± 1.0 
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FIG. 4. Fiducial acceptances mapping each measured bin of the p^, \y^\, A^jets, and pi^, distributions from the inclusive 
phase space to the respective fiducial regions. The factors are derived using Powheg for ggF and VBF production and 
Pythias for VH, tiH, and hbH. The width of the band indicates the uncertainty from five sources: missing higher order 
corrections, PDF variations, changing the Higgs boson mass and the production mode composition, and variations of the 
hadronization/underlying event tunes (see main text for further details). The PDF uncertainty is the largest individual 
contribution to the total uncertainty. 
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FIG. 5. Absolute and fractional cross sections in bins of jet multiplicity for inclusive Higgs boson production at ^/s = 8 TeV 
measured by combining the H > 77 and H —>■ ZZ* —>■ analyses using 20.3 fb“^ of pp collisions. The top plot show the cross 
section in inclusive jet bins, while the other plots have exclusive jet binning (except for the > 3 jets bin). In the lower two plots, 
the cross sections of the two channels are shown individually, defined by the fiducial cross sections corrected for acceptance and 
branching ratio (see Eq. 1). These cross sections have partially correlated systematic uncertainties that are considered in the 
combined measurement. 
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pH [GeV] 



FIG. 6. Differential cross sections (left) and shapes (right) of the Higgs boson transverse momentum (top), absolute rapidity 
(middle) and leading jet transverse momentum (bottom) of inclusive Higgs boson production at ^/s = 8 TeV measured in the 
H > 77 and H —>■ ZZ* —>■ A£ final states using 20.3 fb“^ of pp collisions. Both the combined measurements as well as the 
individual channels are shown. 
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Result tables 

Tables IV-VII present the measured differential cross 
sections and Tables VIII-XI report the corresponding 
shape measurements. 

TABLE IV. Measured cross section in bins of p^. The first 
uncertainty is statistical, the second is systematic. 

Bins [GeV] da/dp^ [pb/GeV] 


0-20 

0.20 ± 

0.15 ± 

0.01 

20-30 

0.88 ± 

0.27 ± 

0.04 

30-40 

0.46 ± 

0.23 ± 

0.02 

40-50 

0.56 ± 

0.20 ± 

0.03 

50-60 

0.28 ± 

0.18 ± 

0.01 

60-80 

0.188 ± 

0.090 ± 

0.009 

80-100 

0.136 ± 

0.067 ± 

0.006 


100-200 0.0214 ± 0.0091 ± 0.0010 


TABLE V. Measured cross section in bins of \y^\. The first 
uncertainty is statistical, the second is systematic. 

Bins da/d\y^\ [pb] 

0.0-0.3 15.3 ± 5.4 ± 0.7 

0.3-0.6 14.5 ± 5.6 ± 0.6 

0.6-0.9 10.8 ± 5.3 ± 0.5 

0.9-1.2 22.3 ±6.7 ±1.0 
1.2-1.6 18.8 ± 7.9 ± 0.9 
1.6-2.4 9.4 ± 5.3 ± 0.4 


TABLE VI. Measured cross section in bins of Ajets. The first 
uncertainty is statistical, the second is systematic. 

Bins da/dA'jeta [pb] 

0 15.3 ± 4.3 ± 0.8 

1 11.9 ± 2.7 ± 0.6 

2 4.1 ± 1.5 ± 0.3 
> 3 2.57 ±0.97 ±0.32 


TABLE VIII. Measured fractions in bins of pij?. The first 
uncertainty is statistical, the second is systematic. 


Bins [GeV] 

l/cr da/dp? [1/GeV] 

0-20 

0.0055 ± 

0.0036 ± 

0.0006 

20-30 

0.0258 ± 

0.0077 ± 

0.0003 

30-40 

0.0133 ± 

0.0068 ± 

0.0002 

40-50 

0.0166 ± 

0.0060 ± 

0.0002 

50-60 

0.0079 ± 

0.0053 ± 

0.0001 

60-80 

0.0055 ± 

0.0027 ± 

0.0001 

80-100 

0.0041 ± 

0.0021 ± 

0.0000 


100-200 0.00060 ± 0.00026 ± 0.00001 


TABLE IX. Measured fractions in bins of The first un¬ 
certainty is statistical, the second is systematic. 

Bins 1/(T d(7/d\y^\ 

0.0-0.3 0.46 ±0.16 ±0.00 

0.3-0.6 0.43 ±0.16 ±0.00 

0.6-0.9 0.32 ± 0.15 ± 0.00 

0.9-1.2 0.66 ±0.19 ±0.00 

1.2-1.6 0.57 ± 0.23 ± 0.00 

1.6-2.4 0.27 ± 0.13 ± 0.00 


TABLE X. Measured fractions in bins of Ajets- The first un¬ 
certainty is statistical, the second is systematic. 

Bins l/cr dcr/dVjets 
0 0.447 ± 0.078 ± 0.010 

1 0.353 ± 0.071 ± 0.005 

2 0.123 ± 0.043 ± 0.003 
> 3 0.077 ± 0.029 ± 0.007 


TABLE VII. Measured cross section in bins of piji. The first 
uncertainty is statistical, the second is systematic. . 

Bins [GeV] da/dpif[ [pb/GeV] 

0-30 0.51 ± 0.14 ± 0.03 
30-50 0.36 ± 0.11 ± 0.02 
50-70 0.156 ± 0.069 ± 0.009 
70-100 0.111 ± 0.050 ± 0.006 
100-140 0.055 ± 0.026 ± 0.004 


TABLE XL Measured fractions in bins of pi^. The hrst un¬ 
certainty is statistical, the second is systematic. 

Bins [GeV] TJ7do]&^^'[l/GeW] 


0-30 

0.0162 ± 

0.0027 ± 

0.0003 

30-50 

0.0117 ± 

0.0032 ± 

0.0002 

50-70 

0.0051 ± 

0.0022 ± 

0.0001 

70-100 

0.0035 ± 

0.0016 ± 

0.0001 


100-140 0.00180 ± 0.00088 ± 0.00007 
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Uncertainty correlation tables 

Tables XII-XV contain the correlation matrices of the 
differential cross section measurements and Tables XVI- 
XIX those of the differential shape measurements. 


TABLE XII. Correlation matrix for the total uncertainty of 
the differential cross-section measurement in bins of p^. 



Bin 1 

Bin 2 

Bin 3 

Bin 4 

Bin 5 

Bin 6 

Bin 7 Bin 8 

Bin 1 

1.00 

0.01 

0.01 

0.01 

0.00 

0.01 

0.01 

0.01 

Bin 2 

0.01 

1.00 

-0.13 

-0.10 

0.01 

0.01 

0.01 

0.02 

Bin 3 

0.01 

-0.13 

1.00 

-0.08 

0.01 

0.01 

0.01 

0.01 

Bin 4 

0.01 

-0.10 

-0.08 

1.00 

0.01 

0.01 

0.01 

0.01 

Bin 5 

0.00 

0.01 

0.01 

0.01 

1.00 

-0.16 

-0.11 

0.01 

Bin 6 

0.01 

0.01 

0.01 

0.01 

-0.16 

1.00 

-0.11 

0.01 

Bin 7 

0.01 

0.01 

0.01 

0.01 

-0.11 

-0.11 

1.00 

0.01 

Bin 8 

0.01 

0.02 

0.01 

0.01 

0.01 

0.01 

0.01 

1.00 


TABLE XIII. Correlation matrix for the total uncertainty of 
the differential cross-section measurement in bins of 



Bin 1 

Bin 2 Bin 3 

Bin 4 

Bin 5 

Bin 6 

Bin 1 

1.00 

0.01 

0.01 

0.02 

0.01 

0.01 

Bin 2 

0.01 

1.00 

0.01 

0.02 

0.01 

0.01 

Bin 3 

0.01 

0.01 

1.00 

0.01 

0.01 

0.01 

Bin 4 

0.02 

0.02 

0.01 

1.00 

0.02 

0.01 

Bin 5 

0.01 

0.01 

0.01 

0.02 

1.00 

-0.28 

Bin 6 

0.01 

0.01 

0.01 

0.01 

-0.28 

1.00 


TABLE XIV. Correlation matrix for the total uncertainty of 
the differential cross-section measurement in bins of Njets- 




Bin 1 

Bin 2 

Bin 3 

Bin 4 

Bin 

1 

1.00 

0.03 

-0.02 

0.02 

Bin 

2 

0.03 

1.00 

-0.04 

0.05 

Bin 

3 

-0.02 

-0.04 

1.00 

-0.04 

Bin 

4 

0.02 

0.05 

-0.04 

1.00 


TABLE XV. Correlation matrix for the total uncertainty of 
the differential cross-section measurement in bins of pijl. 




Bin 1 

Bin 2 

Bin 3 

Bin 4 

Bin 5 

Bin 

1 

1.00 

0.03 

0.02 

0.02 

0.01 

Bin 

2 

0.03 

1.00 

0.02 

0.02 

0.02 

Bin 

3 

0.02 

0.02 

1.00 

0.02 

0.01 

Bin 

4 

0.02 

0.02 

0.02 

1.00 

-0.19 

Bin 

5 

0.01 

0.02 

0.01 

-0.19 

1.00 


TABLE XVI. Correlation matrix for the total uncertainty of 
the differential shape measurement in bins of pJjl. 




Bin 1 

Bin 2 

Bin 3 

Bin 4 

Bin 5 

Bin 6 

Bin 7 

Bin 8 

Bin 

1 

1.00 

-0.28 

-0.22 

-0.18 

-0.16 

-0.17 

-0.11 

-0.12 

Bin 

2 

-0.28 

1.00 

-0.34 

-0.28 

-0.11 

-0.11 

-0.08 

-0.08 

Bin 

3 

-0.22 

-0.34 

1.00 

-0.21 

-0.08 

-0.09 

-0.06 

-0.06 

Bin 

4 

-0.18 

-0.28 

-0.21 

1.00 

-0.07 

-0.07 

-0.05 

-0.05 

Bin 

5 

-0.16 

-0.11 

-0.08 

-0.07 

1.00 

-0.26 

-0.18 

-0.04 

Bin 

6 

-0.17 

-0.11 

-0.09 

-0.07 

-0.26 

1.00 

-0.18 

-0.05 

Bin 

7 

-0.11 

-0.08 

-0.06 

-0.05 

-0.18 

-0.18 

1.00 

-0.03 

Bin 

8 

-0.12 

-0.08 

-0.06 

-0.05 

-0.04 

-0.05 

-0.03 

1.00 


TABLE XVII. Correlation matrix for the total uncertainty of 
the differential shape measurement in bins of |j/^|. 




Bin 1 

Bin 2 

Bin 3 

Bin 4 

Bin 5 

Bin 6 

Bin 

1 

1.00 

-0.09 

-0.09 

-0.11 

-0.10 

-0.21 

Bin 

2 

-0.09 

1.00 

-0.09 

-0.12 

-0.11 

-0.22 

Bin 

3 

-0.09 

-0.09 

1.00 

-0.11 

-0.10 

-0.21 

Bin 

4 

-0.11 

-0.12 

-0.11 

1.00 

-0.13 

-0.27 

Bin 

5 

-0.10 

-0.11 

-0.10 

-0.13 

1.00 

-0.66 

Bin 

6 

-0.21 

-0.22 

-0.21 

-0.27 

-0.66 

1.00 


TABLE XVIII. Correlation matrix for the total uncertainty 
of the differential shape measurement in bins of Ajeta. 




Bin 1 

Bin 2 

Bin 3 

Bin 4 

Bin 

1 

1.00 

-0.77 

-0.37 

-0.26 

Bin 

2 

-0.77 

1.00 

-0.16 

-0.11 

Bin 

3 

-0.37 

-0.16 

1.00 

-0.05 

Bin 

4 

-0.26 

-0.11 

-0.05 

1.00 


TABLE XIX. Correlation matrix for the total uncertainty of 
the differential shape measurement in bins of pi,!. 

Bin 1 Bin 2 Bin 3 Bin 4 Bin 5 

Bin 1 1.00 -0.63 -0.38 -0.36 -0.20 

Bin 2 -0.63 1.00 -0.12 -0.11 -0.06 

Bin 3 -0.38 -0.12 1.00 -0.07 -0.04 

Bin 4 -0.36 -0.11 -0.07 1.00 -0.28 

Bin 5 -0.20 -0.06 -0.04 -0.28 1.00 
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Gluon fusion cross section 

Figure 7 shows the measurement of the Higgs boson 
production cross section compared to a range of theory 
predictions, including LHC-XS, the result used by the 
ATLAS and CMS collaboration in Run 1, for which the 
ggF part is accurate to NNLO+NNLL in QCD [10], as 
well as ggF cross section calculations that attempt to 
go beyond NNLO, including the recently completed full 
N^LO prediction. Details about the various predictions 
are presented in Table XX, and the central values and a 
breakdown of the uncertainties of the calculations as well 
as the measurement are reported in Table XXL 


J2 

I 50 
T 
& 

O 45 
40 
35 
30 
25 


20 
15 

FIG. 7. Measured total cross section of Higgs boson produc¬ 
tion compared to different theoretical calculations. 


ATLAS \li = 8TeV, 20.3 fb' 
^ 5 H^ZZ*—>4/ 

} comb, data ■ syst. one. 


NNLO NPLO NNLO+a-2 NNLO+^ii approx. N“LO approx. N“LO 

+NNLLlhr. +NNLLthr. +N^LLthr. 

Data LHC-XS ADDFGHLM ABNY STWZ dFMMV BBFMR 



125.4 GeV 

XH = VBF +VH+ ftH + bbH 
El QCD scale uncertainty 
E QCD scale and N^LO approx, uncert. 

■ Tot. uncert. (scale, N^LO approx. ® PDF+a^) 


TABLE XX. Summary of the ggF predictions used in the 
comparison with the measured cross sections. The second 
column states the order in QCD perturbation theory and 
which threshold resummation is applied, if any. Further de¬ 
tails are provided in the footnotes. All predictions are for 
niH = 125.4 GeV and ^/s = 8 TeV. 


Total cross- 

-section calculations 

LHC-XS [10] 

NNLO±NNLL^±’‘= 

ADDFGHLM [27-30] N+0‘‘±’‘= 

ABNY [47] 

NNLO±NNLL‘‘±’‘=’‘^’® 

STWZ [31] 

NNLO'^± 

dFMMV [48] 

approx. N^LO"^ 

BBFMR [49-51] 

approx. N+0±N+L‘‘±’‘= 


“ Considers b- (and c-) quark masses in the gg ^ H loop 
Includes electroweak corrections 
Based on MSTW2008nnlo [18] (oa from PDF set) 

Uses TT^-resummed gg ^ H form factor 
® In the counting of Ref. [47], the result has N^LL accuracy 

For the predictions, uncertainties from renormaliza¬ 
tion, factorization and, where appropriate, resumma¬ 
tion scale variations as well as uncertainties due to ap¬ 
proximation or missing terms beyond NNLO are pro- 


TABLE XXL Central values and uncertainties for the differ¬ 
ent ggF predictions and the data. 


Name 


o-gg-s-H [Pb] 


Data-AH"^ 


30.0 ±5.3 (stat) ±1.6 (sys) 


LHC-XS 19.15 +\il (scale) 
ADDFGHLM 20.55 (scale) 


ABNY 

STWZ 

dFMMV 

BBFMR 


+ 1.44 
-1.32 

+ 1.60 
-1.44 


19.54 155 ®4 (scale) 
20.41 ±1.18 (scale) 
21.12 +1-11 (scale) 
21.32 (scale) 


(pdf) 

(pdf) 

^1.35 (pdf) ±0-78 (appr.) 

tni? (pdf) 

il'.ie (pdf) ±0.56 (appr.) 
-\'a 7 (pdf) ±1-39 (appr.) 


^ Non-ggF cross section 

(TxH = 3-Ollto 06 ^ 0.09 (pdf) pb, subtracted from the 

measured inclusive cross section: 33.0 ±5.3 (stat) ± 1.6 (sys) pb. 


vided separately for each prediction. The same rela- 
five PDF uncertainty of assigned to all ggF 

predictions, except for the ADDFGHLM prediction for 
which this uncertainty is increased to correspond¬ 

ing to the change in MSTW2008nnlo uncertainty ob¬ 
served by the group when changing the matrix element 
from NNLO to N^LO. The non-ggF contribution {axH = 
d-Olfofl (scale) ± 0.09 (pdf) pb, XH = VBF + VH + 
ttH ± bbH) is added to the ggF predictions to be able 
to compare to the data in Fig. 7. 

As detailed in Table XX, all inclusive predictions use 
the same PDF set but differ in the perturbative calcula¬ 
tion. Four of the predictions apply both electroweak cor¬ 
rections and consider finite b- and c-quark masses. These 
corrections have non-negligible impacts on the ggF cross 
section; the electroweak correction results in an increase 
of approximately 5%, while the bottom and charm cor¬ 
rections give a 0(5 — 10%) reduction depending on pre¬ 
cisely how they have been implemented in the calcula¬ 
tions. They therefore have an opposite effect on the total 
cross section such that their numerical effects partially 
cancel. The STWZ and dFMMV predictions consider 
neither of these corrections. 

The calculations take different approaches to approx¬ 
imately evaluate the ggF cross section beyond NNLO. 
Therefore the preferred scale for each calculation differs, 
and the choice of scale and the precise scale variations 
applied was left to the authors of the calculations. The 
LHC-XS, ABNY, STWZ, and BBFMR predictions use 
a central scale of no = tuh as their overall scale, while 
dFMMV and ADDFGHLM use = mu 1^- 
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Compatibility between predictions and data 

Tables XXII and XXIII present compatibility tests be¬ 
tween the differential predictions and the measured cross 
sections and shapes respectively. The theory uncertain¬ 
ties are assumed to be Gaussian and to be fully correlated 
between bins. 


TABLE XXII. p-values quantifying the compatibility between 
predictions and the data for the differential cross sections. 
The theory uncertainties are assumed to be Gaussian and to 
be fully correlated between bins. 



P? \y^\ 

pir' 

HRes 

2% 14% 

- 

STWZ 

- 

26% 

JetVHeto 

- 

24% 


TABLE XXIII. p-values quantifying the compatibility be¬ 
tween predictions and data for the differential shapes. The 
theory uncertainties are assumed to be Gaussian distributed 
and fully correlated between bins. 




Pt 


Pt 

HRes 

15% 

64% 

- 

NNLOPS 

10% 

64% 

64% 

SHERPA 2.1.1 

22% 

63% 

88% 

MG5 aMG@NLO 

8% 

60% 

88% 


Non-perturbative correction factors 

Table XXIV presents multiplicative non-perturbative 
correction factors and associated uncertainties that are 
applied to correct analytical parton-level predictions pre¬ 
sented in this Letter to particle level. These corrections 
account for hadronization and multiple parton interac¬ 
tions, and are derived based on a number of underlying 
event and showering tunes applied to the Higgs boson 
production MC samples used in the analysis. 


TABLE XXIV. Non-perturbative factors in percent with sys¬ 
tematic uncertainties, accounting for the impact of hadroniza¬ 
tion and underlying event. 


Bin 


1 


2 



3 


4 


5 


6 

7 

8 

Pt 


99.5 ± 

1.0 

100.5 ± 

1.3 

99.9 

± 

1.1 

99.9 ± 

1.4 

100.2 ± 

1.4 99.9 ± 1.1 

100.4 ± 1.5 

100.0 ± 1.0 

\y^\ 


100.0 ± 

0.8 

100.1 ± 

0.7 

99.8 

± 

0.7 

100.1 ± 

0.8 

100.0 ± 

0.7 99.8 ±0.6 



Vjets 

, excl 

100.6 ± 

2.1 

99.8 ± 

2.9 

99.3 

± 

5.0 

96.7 ± 

6.7 






Vjets 

, incl 

100.0 ± 

0.0 

99.2 ± 

2.5 

98.7 

± 

5.3 

96.8 ± 

6.7 






Pt 


100.6 ± 

2.1 

99.6 ± 

2.6 

99.2 

± 

2.9 

99.1 ± 

3.2 

98.7 ± 

3.6 




















